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In this note, we calculate the branching ratio and CP asymmetry parameters of Bg np in the 
framework of perturbative QCD approach based on kr factorization. This decay can occur only 
via annihilation diagrams in the Standard Model. We find that (a) the charge averaged Br{Bs 
TT+p" +7r"p+) is about (9 - 12) x 10"^; Br{Bs 7r°p°) ~ 4 x 10~^; and (b) there are sizable CP 
asymmetries in the processes, which can be tested in the near future Large Hadron Collider beauty 
experiments (LHC-b) at CERN. 
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The B meson rare decays provide a good place for testing the Standard Model (SM), studying CP violation and 
looking for possible new physics beyond the SM. A lot of theoretical studies have been done in recent years, which 
are strongly supported by the running B factories in KEK and Stanford Linear Accelerator Center (SLAC). Looking 
forward to the future CERN Large Hadron Collider beauty experiments (LHC-b) , a large number of Bg and Be mesons 
can also be produced. So the studies of B^ meson rare decays are necessary in the next a few years. 

In this paper, we study the rare decays Bs np in Perturbative QCD approach (PQCD) p^]. Comparing with 
QCD factorization approach PQCD approach can make a reliable calculation for pure annihilation diagrams in 
■ kx factorization. The endpoint singularity occurred in QCD factorization approach can be cured here by the Sudakov 
T— I ' factor from resummation of double logarithms. 

, In PQCD approach, the decay amplitude can be written as: 

in ■ 

5' Amplitude ~ /d^'fcid^fczd'^fcg Tr[C(i)$s ('fci)$,r(fc2)«'p(fc3)i?(fci, fca, /fcs, *)] e~^^*^. (1) 

^ ■ J 

I In our following calculations, the Wilson coefficient C{t), Sudakov factor Si{t){i = i?s, tt, p) and the non-perturbative 
but universal wave function can be found in the Refs. .S,^, .6^, 7J. The hard part H are channel dependent but 
fortunately perturbative calculable, which will be shown below. 

Like the B^ —>■ tt+tt" decay 'si], the Bg — ^ Trp decays are pure annihilation type rare decays, which are difficult to be 
calculated in method other than PQCD approach. Fig. shows the lowest order Feynman diagrams to be calculated 
in PQCD approach where the big dots denote the quark currents in four quark operators for Bs t:^ p^ according 
^ • to the effective hamiltonian of b quark decay |^. First for the usual factorizable diagram (a) and (b), the sum of 
I (V-A)(V-A) current contributions is given by 

Ma[C] ^ IGTrCFM"^ / dx2dx3 / bidbihdh^ {[x30^(a;2)<?!)p(a;3) + 2r^rp(l + X3)0^(a;2)(^;(a;3) 
Jo Jo 

+ 2r,rp(x3 - l)0^'(a;2)0*(x3)]a,(ti)/i„(x2,X3,52,63)exph5,(ti) - Sp{t\)]C{tl) 
- [X2<t>i{x2)(t)p{xz) + 2r^rp(l + X2)<t>^ {x2)4>"p{x3) + 2r^rp(x2 - I)0j(x2)0p(x3)] 

aSl)ha{x;,, X2,bs, 62) eMSAtl) - Sp{tl)]C{tl)}, (2) 

where r^r = "t-ott/^b — ™^/[^7is(™ti + fnd)], Tp — mpjraB- Cp — 4/3 is the group factor of the S'[/(3)c gauge group, 
is light cone distribution amplitude, which describes the momentum distribution of the meson wave function. The 
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function 



X [e{b2 - b3)MMBV^b3)^H^''\MBV^b2) + (62 ^ &3)] (3) 



comes from the Fom'ier transformation of propagators of virtual quark and gluon in the hard part calculations. St{x) 
results from the threshold resummation of QCD radiative corrections to hard amplitudes ^oi^) — Jo{z) + iYQ{z), 
Jq and Yq are Bessel functions. The sum of (V-A)(V+A) current contributions of diagrams (a) and (b) is —Ma[C]. 

For the non-factorizable annihilation diagrams (c) and (d), all three meson wave functions are involved. The 
(V-A)(V-A) operator's contribution is: 

— 1 f°° 

Mc[C] = 647rCFM| / dxidx2dxz I bidbi 62^62 

v'^Nc Jo Jo 

X {[-X2(l>p{x3)(t)^{x2) + rpr^{x3 - X2)(f>l{x3)(t)^ {X2) - rpr^(x2 + X3)(j)p{x3)(j>^ {X2) 

- r^rp{x2 + X3)(f>*p{x3)(j)^{x2) + r^rp{x3 - X2)(j)l{x3)4>^{x2)] 

C{t\)as{t^)h^P {xi,X2.X3MM) eM~SB{t\) ~ S^{t\) - Sp{tl)] - [~X3^p{x3)dpi{x2) 

+ rpr^{x2 - X3)(j)^p{x3)(t)^{x2) - rpr^{2 + X2 + X3)(/>^(x3)0f (X2) 

+ r^rp(2 - X2- X3)(f>l{x3)(j)l{x2) + r^rp{x2 - X3)(j)l{x3)4>^{x2)] 

C{tl)a,{tl)h''i\xuX2,X3MM)eM-SB{tl) - S^{tl) - Sp{tl)]}. (4) 
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For the penguin operators, there are also {V — A){V + A) type operators, whose contribution is given as 



1 f°° 

GATrCpMg / dxidx2dx3 / bidbi b2db2 (f>B{xi,bi) 
''2Nc Jo Jo 



where 



X { [-X3<j)p{x3)(p^{x2) - rpr^{x3 - 0:2 )(/)p(x3) 0^(2:2) - rpr^{x2 + X3)(l>p{x3)(j)^ {X2) 
- r7vrp{x2 + X3)(f>*p{x3)(l)l{x2) ~ r^rp{x3 - X2)(l)l{x3)4>^{x2)] 

C{tl)as{tl)h^P {XUX2,X3MM) eM-SB{tl) - S^tl) ~ Sp{tl)] ~ [-X2M^3)<l^^i^2) 

- rpr„{x2 - X3)(t>l{x3)<P^ {X2) - rprrr{2 + X2 + X3)4>p{x3)(l)^ {X2) 
+ r^rp{2 -X2- X3)(j)\{x3)(j)^{x2) - r^rp{x2 - X3)(j)p{x3)(l>l{x2)] 

C(t2)a,(i2)ft(2)(^^^2.2^2:3,5i,52)exp[-5s(te)-5^(^c)-^p(ic)]}, (5) 



h'^j\xi,X2,X3,bi,b2) = 

0{b2 - bi)^ul^\MBV^b2)MMBV^bi) 



+ (&! 



Ko(AfBF(,)6i 



for Fg.)>0^ 



f hW(Mb 61), for < I ' 



Kq is modified Bessel function and 's are defined by 

F^i) = X1X2 - X2X3; F^2) = xi + X2 + X3 - X1X2 - X2X3. 



(6) 



(7) 



The hard scale t'^s in Eqs.|j21ElEl are chosen as the largest energy scale appearing in each diagram to kill the large 
logarithmic corrections: 



ti 



max(MBV^7 1/^2, 1/&3), 
max(MBy/x^, f 7^2, 1/&3), 

max(AfB^|F2)|, Mb^^, l/bi, f 762) 



max{MB J \F^2)\' ^'^By/x^, f/62 



The total decay amplitude is then 

AiB, ^ TT+p-) = fBMa 

+M, 

and the decay width is expressed as 



(8) 



p 



TT^ P 



(9) 



The most important contribution here is the factorizable penguin diagram, which is CKM enhanced. If we exchange 
the TT and p in Fig. f , by the same method, we can compute the Bs — *■ Tr~ p'^ decay. The expressions are similar. The 
decay amplitude for Bg n^p^ is 

A{Bs ^ = A{Bs ^ n+p-) + AiB, ^ 7r-p+), 

and the decay width can be written as 



(10) 



OlzTT ' 



(11) 



4 



In the following, we first give the branching ratios of Bg np. Just as in Ref. we leave the CKM angle 7 as a free 



parameter in our numerical calculations. Because there are four decay channels: Bg/ Bg — > Tr+p , Bg/Bg — > tt 
is not possible to distinguish the initial state by detecting the final states. We average the sum of Bg/Bg tt^ p" 



, it 

as 



one channel, and Bg/Bg 
refs. 0,13, we get 



TT p"*" as another, which is distinguishable by experiments. Using the same parameters as 



Br{Bg/Bg 
Br{Bg/Bg 
Br{Bg/Bg 



n+p-) = (S.ll^J) X 10-^ 
n-p+) = {bAtlD X 10-^ 
.vrV) = (4.2tr7)xlO-^ 



(12) 



where all channels are averaged for Bg and Bg. 

In Rcf. 11], Beneke etal have estimated the branching ratio of Bg 7r+/5~ in the QCD factorization approach. 
Weak annihilation diagrams are power suppressed in the heavy quark limit and, in general, not calculable in QCD 
factorization approach. In order to avoid the end-poind singularities, they introduced phenomenological parameters 
to replace the divergent integral. With those parameters they estimated that the branching ratio of Bg 7r+p~ 
is (0.03 — 0.14) X 10~^. In PQCD approach, the annihilation amplitude is calculable. In the rest frame of the Bg 
meson, the d or d quark included in p or tt has momentum C'(Mb/4), and the gluon producing them has momentum 
g2 = 0(M^/4). This is a hard gluon, the PQCD can be safely used because of asymptotic freedom of QCD We 
have tested that the bulk of the result comes from the region with ag/w < 0.2, where a figure was show in Ref. [l3l |. 
Our predicted result is larger than their estimation, which can be tested by the future experiments. 



Using the same definition in Ref. [8j, we study the CP violation parameters A^p and 



Bg TT p^ decay. We find the direct CP violation parameter A'Qp{Bg 



in the process of 



is about 2% when 7 is near 100°, 



the small direct CP asymmetry is also a result of small tree level contribution. The mixing CP violation parameter 



{Bg 



ttV) 



is large, whose peak is close to 16% (see Fig.|2l. 
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The CP violations of Bg/B, 
expressed as: 



FIG. 2: Mixing CP violation parameter of Bs yr^p" decay as a function of CKM angle 7 

TT^p^ are very complicated. There are four decay amplitudes, which can be 



h^{7r+p-\H.jf\Bg); 
g^{7r-p+\H,ff\Bg); h ^ p+\H,f f\Bg) . (13) 

We introduce four parameters to describe the CP asymmetries in the processes, which are given by 0| 



g={iT+p-\H,ff\Bg)- 



C-rrfi 



1/ 



1 



£ + £' )l 



(14) 
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where 

_ \g?-\h\^ _ -2Im{h/g) 

\h?-\9? -2Im{g/h) 
'~\h? + \-9?' '^'^ l + lffAP ■ ^ ^ 

We calculate the above four CP asymmetry parameters and show the CKM angle 7 dependence in Fig. 13 We do 
not plot the C-r^p in this figure, since its value is near zero. The decay branching ratios depends heavily on the shape 
of wave functions and decay constants etc. But the CP asymmetry should not since the dependence will be cancelled. 
The direct CP asymmetry can be affected by the power corrections and next-to-leading order contributions easily. We 
investigate the CP asymmetry parameters's dependence on the hard scale t in Eq. J^l, which characterize the size of 
next-to- leading order contribution. The CP asymmetry numbers are shown in Table with those uncertainties. By 
changing the hard scale t from 0.9t to 1.3t, we find the CP asymmetries of Bs{Bs) n^p^ change little: for 7 = 60°, 
the uncertainty is less than 1% for C^rp, 4% for AC^rp, 3% for S-^p and 7% for AS'^rp, respectively. The reason is that 
mixing induced CP is dominant here, but the direct CP of Bg — s- tt^/o" really changed much, as shown in Tabled 



TABLE I: CP asymmetry parameters using 7 = 60° with uncertainties 



Scale 






Snp 


^S-TTp 






0.9 t 


0.1% 


73% 


10% 


21% 


0.6% 


11% 


t 


1% 


77% 


11% 


14% 


1.6% 


13% 


1.3 t 


1.8% 


79% 


14% 


10% 


2.6% 


16% 
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FIG. 3: CP violation parameters of B^{Bs) ir^p^ decays: AC (solid line), S (dashed line) and AS (dotted line) as a 
function of CKM angle 7 

In conclusion, we study the branching ratio and CP asymmetries of Bg — > np decays in PQCD approach. Wc find 
the branching ratio of Bs 7r+p~ +tt~p^ is at order 10~^, which is larger than QCD factorization's estimation. We 
also predict CP asymmetries in the process, which may be measured in the future LHC-b experiments. 
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